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Abstract 



' We study on the lattice the correlator of heavy-quark currents in the vicinity of vanishing 



\ momentum. The renormalized charmed quark mass, the renormalized strong coupling 

constant and gluon condensate can be defined in terms of the derivatives of that correla- 
tor at zero momentum. We analyze quenched Monte-Carlo data on a small lattice 8^ * 16 
for /5 = 6. We generalize dispersion relations to the lattice theory in a simple way and 
use them successfully to fit the correlator at both small and large distances. We fit the 
short- distance part of the correlator with the relevant expressions of perturbative QCD on 
the lattice and obtain the value of the renormalized quark mass m^^^{m^ = 1.20(4) GeV. 
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Introduction 

At present the accurate calculations of the strong coupling constant Jll and the masses 
of heavy quarks , are not only the subject of intensive studies in the continuum QCD. 
Monte-Carlo Lattice simulations seem to provide us with competitive ways to evaluate 
as and heavy-quark masses P] with considerable accuracy. A systematic study 

of the correlator of heavy-quark currents n(g^) near zero momentum = on the 
lattice was suggested recently as a way to calculate the renormalized parameters of 
perturbative QCD, the strong coupling constant as and heavy quark mass irtc-, as well as 
the main parameter of non-perturbative QCD, the vacuum gluon condensate C-^-* = < 



In this paper we formulate the method in more detailes and analyze Monte-Carlo data 
on relatively small lattice 8^*16, as a pi/o^ study. We concentrate on the renormalized mass 
of the charm quark rric. The renormalized parameters of the QCD Lagrangian, 
and mc(/i^) serve as the boundary conditions to the renormalization group equations. 
Here /x^ is some physical normalization point, which stays finite as the cut-off is removed, 
say the lattice spacing a — > 0. This choice of a physical normalization point allows one 
to easily make use of phenomenological data in the determination of Aqcd and current 
quark masses on the Lattice. It also helps to isolate such non-perturbative quantities as 
the gluon condensate, whose naive definition would imply a strong power-like dependence 
on the normalization point This is the main difference between the proposal 



and the works ||T2| where the computation of the gluon condensate on the lattice was 
pioneered. 

Because of the asymptotic freedom the boundary conditions to the renormalization 
group equations in QCD should be imposed at large Euclidean momenta corresponding 
to small space-like distances. The renormalized parameters of QCD are essentially the 



short- distance quantities. We would like to emphasize the difference between the pole- 
mass of quarks and Euclidean quark mass. The pole mass is a parameter of the on-shell- 
subtraction scheme: the ultraviolet singularities are subtracted on the quark mass-shell. 
This scheme does not appeal because of the complicated infrared properties of QCD: 
quarks do not exist on the mass-shell. The pole mass of quarks is well defined only within 
the first orders of perturbation theory, while the Euclidean mass is defined in the domain 
where perturbative QCD apphes. In the popular these days language one would say that 
the pole-mass suffers from the infrared-renormalon ambiguities while the Euclidean mass 
does not. 

One way to calculate the Euclidean mass on the lattice was considered in 0. The 
prescription is: a) within the lattice theory calculate the pole mass m* of the quark 
propagator as a function of the parameters of the lattice theory m*[/t, r, /3, ...] pertur- 
batively in powers of the inverse bare coupling constant (3. Here : (3 = Q/g"^ and K,r 
are parameters of the Wilson propagator. This perturbative expansion is good at large 
/3, where the asymptotic scaling holds: the function g'^{a) solves the perturbative renor- 
malization group equation. Volume has been always considered infinite in these analytic 
calculations. The essential ingredient of this calculation is the bare quark mass m(a) , 
which is obtained by fitting the Monte-Carlo data for the pion mass with the ansatz: 
m(a) = — l//€cr)/2. The motivation for this fit is that ~ triquark is expected in 

the chiral limt. The fit is done for ks corresponding to light fermions in the finite volume 
and for relatively small /5s, where (3 does not scale, h) Use the relation of the continuum 
theory between the quark pole mass and the Euclidean mass m*[m(/i), ^^(yu)]. 

Our approach is somewhat alternative to . We use the correlator of hadronic cur- 
rents of heavy quarks to determine the Euclidean quark mass. We never refer to the pole 
mass. The correlator of the heavy-quark currents is reliably calculated in perturbative 
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QCD at vanishing momentum, because that point is away from the nearest threshold due 
to heavy quark-antiquark pair production. The heavy-quark threshold 4m^ is the relevant 
virtuality here. The perturbative loop-expansion is built in powers of small parameter 
as(4m^). One loop of free quarks with the renormalized Euclidean mass gives a very good 
approximation to the correlator of the charm-quark currents near vanishing momentum 



| T0|| . This is in contrast to the correlator of light quarks, which is not analytic at vanishing 
momentum and cannot be computed at that point within perturbative QCD. 

The correlator of the gauge-invariant local hadronic currents 11 (x) is a very conven- 
tional thing to study in Lattice QCD. The long-distance behavior of the correlator is 
determined by the low-lying resonance. To obtain the resonance mass one fits the zero- 
momentum correlator n(t) by cosh[{t — L / 2)mres\i which describes propagation of a single 
particle of mass irises in a finite volume L in the continuum theory (a = 0). We do the 
same thing: fit the long-distance part of the correlator of heavy-quark currents with the 
single-resonance approximation to that correlator. This allows us to fix the scale of our 
Monte-Carlo data. In addition to that we study the short-distance part of the correlator 
as well. We fit it with the corresponding expressions of perturbative QCD in order to 
extract the fundamental parameters of perturbative QCD: the renormalized quark mass 
and the strong coupling constant. 

In order to separate the short-distance part of the correlator from its long-distance part 
we study derivatives of the correlator n(g^) with respect to momentum at the origin 
= 0, called moments [jlOl- The lower derivatives probe the vicinity of the origin, which 
is 4m^ GeV^ away from the charm quark-antiquark threshold, hence it is in the domain of 
aymptotic freedom. The lower moments are reliably computed in perturbative QCD. We 
use Monte-Carlo data for the lower moments to extract the parameters of perturbative 
QCD. Higher moments probe the values of away from the origin. They represent the 
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long-distance part of the correlator and are entirely saturated by the low-lying resonance. 
We use those moments in order to fix the scale of the lattice theory. 

The next section 1.1 is a short introduction into the so-called dispersion theory of 



charmonium |TT[. We describe analytic properties of the two-point correlator of charm- 
quark currents in the continuum theory and explain why it is reliably calculated near 
vanishing momentum in perturbative QCD. The other sections are our attempt to extend 
the dispersion theory of charmonium to the Lattice QCD. 

In section 1.2 we discuss an extension of the dispersion relations to the lattice theory. 
Since we deal with the lattice theory of heavy quarks significant cutoff effects are expected. 
At /? = 6 the lattice spacing a ~ (2 GeV)~^. For the charm quark l.SGeV one 

has the parameter of the deviation of the fermionic action from the continuum a-rric ~ 0.6 
- moderately, but not comfortably small. Actually one would like to probe virtualities 
/i^ ~ 4ml lattice, then a- fi 0{1). From the phenomenological point of view we 

are trying to generate a resonance of mass ~ ttlj/^ ~ S.lGeV. Then the resonance mass 
in the units of the lattice spacing a ■ nires ~ 1-5. Such lattice theory is obviously coarse 
enough to exhibite certain cutoff effects in the propagation of a single resonance. We 
explore the cutoff effects via fitting Monte-Carlo data with the following generalization of 
the dispersion relation for the two-point correlator Il{q^) to the lattice theory: 

fds^ =^ Us — , . (1) 

J s + J s + ^E>m'(gW2) 

where L represents the finite volume of the box, and the function p{s, a, L) approaches the 
corresponding spectral density of the continuum theory in the limit a — 0, L — > oo. We 
find significant evidence that although the l.h.s. of (|l]) exhibites strong deviations due to 
a 7^ & L <^ cxD, the dependence of the spectral density on a and L is rather weak on the 
lattice under consideration. In other words, most of the cutoff effects in the heavy-quark 
systems are successfully taken into account by the naive modification of the denominator 



in i). 

In section 2 we compare the moments of the correlator in a single-resonance approxi- 
mation on the finite lattice and in the continuum theory. We compare the fit based on the 
dispersion relation (|l]) with the traditional fit to cosh[{x — L/2)mres] for different values 
of the parameter a ■ rrires- We analyze Monte-Carlo data obtained on a blocked lattice 
8^ X 16 with ^ 500MeV, which is supposed to show correctly the long-distance part 
of the correlator, where the single resonance dominates. We find perfect fit to the data 
with help of eqn.(|lD using one fit-parameter, the resonance mass: p(s) = 6{s — m^.^^). 

In order to see the manifestation of the short-distance phenomena, we have generated 
Monte-Carlo configurations on a small (unblocked) 8^ x 16 lattice at /3 = 6, corresponding 
to ~ 2GeV. We use quenched propagators of the clover- and tadpole-improved Wilson 
fermionic action. We find that the long-distance part of the correlator (higher moments) is 
still well fitted with a single resonance. We reproduce the masses of the pseudoscalar (r/c), 
vector (J/ip) and scalar (xco) mesons of the charmonium spectrum with quite reasonable 
accuracy. For lower moments (~ short-distances) though, Monte-Carlo data deviate from 
the single-resonance curve. However, that deviation proves to be described very well by 
the contribution of the excited states - a smooth hadronic continuum spectrum, which we 
incorporate via the dispersion relation ([l|). 

In section 3 we study the lower moments of the correlator, corresponding to the short- 
distance physics. They must be reliably reproduced in perturbative QCD, which means 
that one-loop of free quarks gives the dominant contribution. The mass of those quarks 
is the renormalized heavy-quark mass by construction. Since there is a reliable loop 
expansion for the lower moments of the correlator of heavy-quark currents one can study 
systematic deviations of the coefficients of that expansion order by order. 

First, we fit the Monte-Carlo data for the lower moments by the correlator composed 
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of free Wilson quark propagators. Those propagators correspond to the quarks of mass 
rric = 1.37GeV in the continuum theory. To be more accurate we study a relation 
between the mass of the free Wilson action and the corresponding fermion mass of the 
continuum theory in the finite volume because the Monte-Carlo is done in the finite (and 
relatively small) volume. Again, we find evidence that the dispersion relation (|l]) works 
well, this time - for the continuum spectrum of quark-antiquark pair. To specify the 
subtraction scheme of the renormalized quark mass we take into account the two-loop 
ag-correction. With this refined treatment we obtain in the minimal subtraction scheme 
m^^'^(mc) = 1.22{2)GeV from the vector channel and m^^{m^ = 1.18{2)GeV from the 
pseudoscalar channel. We attribute the difference in these two values to the absence of 
scaling on the given lattice. 

1 Correlator of the heavy-quark currents 

1.1 Dispersion theory of charmonium 

Consider the correlator of vector currents of charmed quarks: 

n(g')^. = iJdxe''>^<0\T{j^{x)UO)}\0> (2) 

where = cT c. The correlator (|^) in the continuum theory satisfies the standard 
dispersion relation: 

/p(s^ 

^'^ g2(g ^ g2) + + d2iq^qu - g^^^q^) (3) 

The correlator n(g^)^^,y is defined up to the polynomial, which coefficients di^2 have ultra- 
violet singularities. The polynomial with the coefficients di and d2 is a local part of the 
correlator. In the X-space one has: 

^{^tf = d,g,,,6^'\x) + d2 {d,d, - g,.d') 6^'\x) (4) 
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is essentially electromagnetic current of charmed quarks. The correlator de- 
scribes electromagnetic vacuum polarization due to strong interactions. The coefficient 
di renormalizes photon mass, while the coefficient ^2 renormalizes photon wave function. 
If the ultraviolet singularities are subtracted on the photon mass shell = one has 
di = d2 = for the renormalized correlator. 

The spectral density p(s) is proportional to the inclusive cross-section of the produc- 
tion of charmed and anticharmed hadrons in e"'"e~-annihilation via the electromagnetic 
current of heavy quarks. From experiment one knows that p{s) has a large peak due to 
J/ip-meson. Higher radial excitations with smaller contribution to the spectral density 
overlap with the continuum spectrum of multiparticle states. The simple and phenomeno- 
logically very successful ansatz for p(s): 

contains a single low-lying resonance of mass m^es = ^T^j/ip-, residue /, proportional to the 
electrwomagnetic width of the J/^/^-meson and a smooth continuum spectrum with some 
effective threshold So > fnf.^^. The coefficient l/An'^ comes from one loop of free quarks. 
The fact that continuum hadronic spectrum at high energies s > is identical to the 
perturbative spectrum of free quarks (with small radiative corrections) is usually referred 
to as global duality This property is a consequence of asymptotic freedom. 

Since the poles and discontinuity of the correlator Il{q'^)^i, are located away from the 
origin = the correlator is analytic at = 0. Following define moments of the 
polarization operator (0) as: 

If/ d ' 



■Mn = -{{-—] n,,(g^) (6) 

) g2=0 



n\ [\ dqy 

From the spectral representation one obtains: 



8 



Obviously the first two moments A^o,i suffer from ultraviolet singularities wheras all the 
other moments are physical quantities: the ultraviolet singularities in them can be ob- 
sorbed via the renormalization of the Lagrangian parameters, the strong coupling constant 
and the charmed quark mass. 

The quantity we study below is the ratio r„ of the neighboring moments: r„ = 
M-n+i/ M-n ■ The phenomenological spectrum (|^) has the following ratios 

1 '^n .ON 



^les Vn-1 



2 \ " 



One can see that the continuum contribution to the higher moments rapidly decreases. 
Higher moments are saturated by large distances, where the single low-lying resonance 
dominates. The ratios (|^) corresponding to the single-resonance approximation to the 
correlator (0) are particularly simple: 

rT = (9) 

where m^es stands for the resonance mass. This is illustrated in Fig.l. 

Lower moments receive substantial contribution from the continuum spectrum. They 
probe the vicinity of the origin = 0. From the QCD point of view that region is in the 
domain of asymptotic freedom: the nearest threshold due to a pair of charmed quark and 
antiquark is away by 4m^ ^ ^qcd- 

The function n(g^) is reliably calculable in perturbative QCD near = 0. Low- 
energy thresholds due to multi-gluon intermediate states, which appear on the 4-loop 
level, are known to give small contributions in the case of heavy quarks [l^ . They may 
be eliminated altogether by giving different masses to the two quarks of the current. 
n(g^ = 0) is a function of the renormalized heavy quark mass mc{fi) and as(/i). The 
relevant renormalization point here, = 2mc, is high enough to ensure the validity of 



the QCD perturbation theory: a;s(4m^) ~ 0.2. Then the apphcabihty of perturbative 
QCD near = imphes the following expansion for the ratios r„: 

an + &„as(4mj - c„ „ (10) 



4m^ \ (Ami 



where {a„, 6„,c„} are known numbers [|TOl- The term ~ a„ comes from one loop of free 
charmed quarks. The term ~ 6„ comes from the two-loop diagrams corresponding to 
one-gluon exchange. The gluon condensate G*^^) = < 0\{as/Tc) G^^^G^j^lO > is the vacuum 
average of the first nontrivial operator in the Wilson operator-product expansion for 
the correlator of heavy-quark currents. The coefficient c„ originating from the Wilson 
coefficient function starts with one loop of the heavy-quark propagators. 

The lower ratios ^2,3,4, originating from short distances, are well reproduced by QCD 
perturbation theory. The typical virtuality corresponding to the nth moment is ~ 4m^/n. 
The coefficient 6„ grows with n. At large n > 8 the perturbation-theory based expan- 
sion (|l^) is irrelevant. The gluon condensate shows up in the intermediate ratios ^"5 6,7 
T0[| , sensitive to larger distances and, hence, nonperturbative fluctuations. It represents 



the first term of the operator product expansion. The operator-product expansion thus 
describes a crossover from the short-distance region to the large- distance region, in the 
vicinity of the domain of asymptotic freedom. The contribution of the gluon condensate 
into the physical quantity r„ is renormalization-scale independent. To extract the value 
of the gluon condensate from r„ one has to know the parameters of perturbative QCD. 
The better one knows and rric from the lower moments ^2,3,4 - the higher the accuracy 
that can be reached in the calculation of the gluon condensate from r^^Qj. 

The phenomenological parameters of hadronic spectrum (|^), based on experimental 
data on the inclusive cross-section of charm-anticharm production in e"'"e~-annihilation. 



are : 



47rV - 0.6, rures ^ 3.1 Gel/, ~ {AGeVf (11) 
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The moment ratios of the corresponding correlator are shown on Fig. 1. One can see that 
the higher moments n > 7 of the complete correlator, shown by bursts, concide with the 
single-resonance approximation to that correlator - the contribution of the J/ip-meson, 
shown by squares. The lower moments, corresponding to the short- distance part of the 
correlator, deviate from the single-resonance line. They are well reproduced with one loop 
of free quarks of mass = 1.26GeV, which moment ratios are shown by diamonds. This 



way the mass of charmed quark was first estimated within QCD in |jTO|. 

In order to make this estimate more accurate one has to know the contribution of con- 
tinuum hadronic spectrum very well, since it is large in the lower moments. However, this 
is demanding with respect to experimental data on the inclusive cross-section of charm- 
anticharm production in the e"'"e~-annihilation. Additional uncertainty here is due to the 
fact that at high enrgies one has to separate the bigger part of that inclusive cross-section, 
coming from the electromagnetic current of heavy quarks, from the smaller part of that 
cross-section, coming from the electromagnetic current of light quarks. In order to avoid 
the uncertainty associated with the hadronic continuum alternative ways to determine 
charmed quark mass have been developed. They are based on the nonrelativistic approx- 
imation to QCD. Potential models are particularly successful for 6-quarks. Charmed 
quarks are light enough to exhibit substantial relativistic effects. In this work we choose 
to stay within the original QCD rather then to study nonrelativistic approximaitons to it. 
In order to avoid the uncertainty due to hadronic contiuum we evaluate the correlator (|^) 
by Monte-Carlo in lattice theory. The lattice calculation requires much less experimental 
data. One only has to know the mass of the low-lying resonance in order to fix the scale 
of the lattice theory. 

The fact that the lower moments are well reproduced in perturbative QCD implies 
an integral relation, called sum rule, between the parameters of perturbative QCD and 
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phenomenological parameters of hadronic spectrum: 



s 



^ PQCnis) = f pphenis) (12) 



where n = {2,3,4} and pqcois) is the spectal density corresponding to the multiloop 
expansion of perturbative QCD (|l^). The parameters of hadronic spectrum, consistent 
with the sum rule ([l2|) and experimental data (available in the case of vector channel), 
are those in eqn.(0). 

One important implication of the sum rules (0) is that smooth hadronic continuum 
spectrum, which form coincides with the corresponding perturbative expression, must be 
incorporated in the model of the hadronic spectral density. At high energies the contin- 
uum spectrum dominates over the contribution of few excited states. The high-energy 
asymptotics of the hadronic spectrum must be the same as predicted by the perturbative 
QCD. Therefore, a finite number of narrow resonances is ruled out as a model of hadronic 
spectral density of the correlator (H). 

One should emphasize that light quarks do not participate in the derivation of the 
spectrum (0) in a sense that the expansion ([T0| ) does not contain diagrams with light 



quarks. Incorporation of light quarks would bring more accuracy, in general, but the 
contribution of light quarks into the lower moments is essentially small, as they appear 
in higher orders of the reliable perturbative expansion in powers of as(4m^ 



1.2 Dispersion relation on the lattice 

In the previous section we have pointed out that the short-distance part of the correlator 
(^"2,3,4) receives substantial contributions from the high-energy (continuum) part of the 
spectrum. When analysing the data on the correlator in lattice QCD we would like to 
know how the correlator with a given spectrum looks like in the lattice theory. How the 
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finite volume and finite cutoff modify the function 11 (x) witli the spectral density as in 
(i? 

We explore the following naive extension of the dispersion relation (|^) to the lattice 
theory: 

n(g^,a,L) = f ds — ■ — 77^ (13) 

In principle, one should expect the spectral density in (0) to have lattice artifacts. In 
practice we find that the simple modification of the denominator as in ( ]T3| ) accounts for 
the major part of systematic deviations due to finite L and a. Eqn.(|13|) allows one to 
fit Monte-Carlo data by a complicated spectrum with several resonances and a smooth 
continuum. 

In the case of a single resonance p{s) = 6{s — m^es) 6qn.(|13|) becomes a familiar 
expression for the lattice propagator of a free field: 



1 



where we have introduced a four- vector n of unit lengh: = qn^, so that Z]^=i nl = 1. 
The extra dependence of the correlator on n,^ is an (9(a^)-effect. 

To see how well eqn.(]T3[) works for the continuum spectrum we fit in section 3 the 
lattice correlator in the free Wilson quarks approximation (which has contuum spectrum) 



to the eqn.(p^) with the quark-antiquark spectral density of the continuum theory. It 
proves to work well on finite lattices of different size. 

Since we are interested in the short- distance behavior of the correlator we have to 
address the question of subtractions. At short distances the correlator is singular. In 
the continuum theory the ultraviolet singularities are proportional to the polynomial of 
finite order in the momentum space. The finite order of that polynomial is a consequence 
of the renormalizability of the theory. The renormalized correlator is the one with the 
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infinite parts of tlie local terms subtracted. In tlie continuum tlieory tlie difference be- 
tween subtracted and unsubtracted correlator is polynomaial in q^. For example, in the 
correlator of the vector currents @) that polynomial is quadratic: di + d2q^. The factor 
here is essentially the inverse Laplacian of the continuum theory. The subtraction 
polynomial contributes only to the first two moments. The other moments do not receive 
any contribution at all from the subtraction polynomial. 

Let us see what the difference is between the subtracted and unsubtracted lattice 
dispersion relation (|T3|). The subtraction: 

1 1 ^ _ ^E>in^(gW2) 

generates non-polynomial terms qf^^ = ^ siii^(?At'2/2) instead of in the numerator. 
To perfom a subtraction in the lattice theory one would have to decompose the correlator 
(|T^) in the following form: 

_ 4 4 

n(g2) = + d2^^sin2(g^a/2) 

+ (^i:sm'iq,a/2)] Ids . (16) 

V« V / s2 (s + ^E;,sm (g^a/2)) 

The factor q^^^ is the inverse of the standard lattice Laplacian: 

/ W) '^'^ ^ ^ ^i^'(^M«/2) = ( 25(^) - ^(^ -a)-5{x + a))/ a" (17) 

Therefore the coefficient d2 in (|16D is precisely the renormalization of the photon wave 
function in the case of vector currents of charmed quarks. The renormalized correlator 



would be given by eqn. (|1^) with rfi = = 0, considered as a function of the renormalized 
Lagrangian parameters. 

The deviation of the gf^j-term from a polynomial q^ is an (9(a^)-cutoff effect. It 
generates 0{a^) systematic deviations in the moments, which decrease rather fast with 
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n. For a single resonance, for example, in the 1- dimensional case one can find analytically 
explicitly: 

In the next sections we will see the size of the contribution of the subtraction term gf^^ 
into the moments on finite lattices. The ratios r„ for n > 3 prove to be practically not 
affected by subtractions. 

2 Moments on a finite lattice 
2.1 Single resonance 

The definition of the moments is designed in such a way that the ratios r„ of the 
neighboring moments of the correlator, saturated by a single resonance, are all equal to 
the inverse mass squared of that resonance. This suggests a simple way to determine 
the resonance mass numerically from given data on the correlator (^. The property (j^) 
however, holds in the continuum theory, that is in the limit a — and L ^ oo. One has 
to see how much the moments of a single resonance are distorted by the finite L and a. 

The extension of the definition of moments (P) to the case of a finite lattice (with 
broken Lorentz invariance) was suggested in ||^: 

■Mn = ,/ , / rfx^x2"n(x) (19) 



22"n!(n+ 1)! 

In the continuum theory eqns. (^ and ([19|) are identical. An alternative to ([T9| ) would be 



to project the correlator (|^) to g* = and define the moments as derivatives with respect 
to q°. In this work we explore the Lorentz invariant definition (p!9|). 

We have evaluated the moments ([1^) of a single resonance (|1^) numerically on a given 



lattice. Figs. 2,3,4 show the moments of a single resonance for different volumes. One can 
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see very strong finite volume effects. For small values of L and arrives the moments are far 
away from the straight line (§) of the continuum theory, illustrated by squares in Fig.l. As 
one increases the volume, the moments approach a straight line corresponding to (P). This 
works better for larger values of arrires (Fig-4), corresponding to more compact resonance- 
objects. The constant behavior of r„ is clearly seen for m^es = 1-5 in the whole range of 
the ratios of interest ri_^7 on a large lattice 32^. For larger n the moments deviate from the 
prediction of the continuum theory @ because higher moments probe larger distances, 
hence they are more sensitive to the boundary of the box. Figure 5 shows the behavior of 
the moments for different resonance masses rUres = 0.5 -r 2 on the lattice 16'^ * 32. Again, 
one can see that finite-volume effects disappear for larger values of the resonance mass in 
the units of the lattice spacing. It is important to realize that these finite-size effects do 
not represent any loss of information. The curves are universal and may (will) be used to 
fit Monte-Carlo data on the same lattices. 

Figs. 2,3,4 show the moments of the subtracted correlator as well. The subtracted 
correlator has vanishing first two moments Aio,i = 0. The ratio ri = Ai2/M.i of 
the subtracted correlator is not to be considered. The two solid lines on Figs. 2, 3, 4, 
corresponding to the subtracted and unsubtracted correlators, are indistinguishable except 
for the second ratio r2. The line corresponding to the subtracted correlator is always the 
one that terminates at n = 2. In the continuum theory the ratio r2 is not affected by 
subtractions because the subtraction polynomial is 0{q^). On the lattice the subtraction 
term qf^^ gives contributions to higher moments. Those contributions are noticable in 
the ratio r2 and negligable in higher ratios. The deviation between the subtracted and 
unsubtracted curve in the ratio r2 is order O(a^) effect. One can see that the deviation 
in the ratio r2 decreases with the resonance mass, because the decrease in the resonance 
mass here essentially implies the decrease in the lattice spacing a. In the limit a —>■ one 
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has gf^j — s> g^, hence the ratio r2 becomes insensitive to subtractions. 

Concentrating on the cutoff effects we extract the resonance mass from the ratio 
r* = rn{L oo), observed on the biggest available lattice 32^, by looking for a plateau 
like the one on the upper curve of Fig. 4 for n < 8. Increasing the volume would extend 
this plateau to larger n, while its position would not change. This would be our procedure 
to fit the resonance mass m^es = (r*)^^/^ via the analysis of the moments. However, even 
if the lattice is large enough to exhibit a perfect plateau for r„ = r*, corresponding to the 
continuum theory, the position of that plateau gives only approximately because 

of finite-a effects, as one can notice by inspecting the upper lines of Figs. 3 and 4. The 
resonance masses m^es are plotted on Fig. 6 versus the input mass m^es- The fit of that 
data gives the following 0(a^)-systematic deviation^: 

1 1 + 0.04 



res 



771^ 771 ^ 

"^res '^res 



(20) 



The fluctuations of the data Fig. 6 around the fit (pOD near rrires = 1 are caused by our 
uncertainty in determining the plateau r„ = r*, which is not well-established yet on a 
32^-lattice for those values of rrires- The value obviously too large to be fitted 

by the (9(a^)-law (|20|) . The systematic correction (^) is rather small for rrires < 1? so 
the present Monte-Carlo data for light quarks (where the resonance masses in units of the 
lattice spacing are typically less than 1) does not see it. 

In Fig. 7 we compare the resonance mass m^es with the conventional way of extracting 
the resonance mass from the correlator (|]) - a fit to cosh{t — L/2) at vanishing three- 
dimensional momentum q = 0. One can see that the resonance mass rfires, extracted from 
the moments is always a better approximation to the expected result rrires- 



^ All the dimensional quantities in the lattice theory are in units of the lattice spacing a throughout 
this paper, although we show the factor a explicitly sometimes to draw attention to the dependence on 
it. 
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2.2 Monte- Carlo data 

In this section we analyze the moments of the correlator (|^), computed by Monte-Carlo in 
the quenched approximation, first on a large but coarse lattice obtained by blocking, then 
on a small but fine lattice. These two studies allow us to check separately our numerical 
control of large and small moments. 

The fermionic action we use Sp = Sw + Ssw has the Wilson term Sw and the clover 



term Ssw 1 19 



Sw = J2{ 



X) 



X 



— I^W 



[^(a;)(r-7^)f/^(x)^(x + /i) + ^(x + /i)(r + 7^)f/;(a;)^(x) ] } (21) 



Ssw = -'^-!^c i){x) F^y{x)a^ui'{,x) (22) 

x,^,v 

where F^y is a field-strengh tensor as the following sum of four plaquettes in the plane 
around point x: 

F,A^) = \il [UuA^) - Uh,A^)] (23) 

The action Sp does not have terms of order 0(a) in the expansion in powers a as a — >■ 
for c = 1 [|l^]. The Green functions computed with the help of action Sf{c = 1) do not 
differ from their limit a — in the order 0{a) provided the fermionic fields ip{x) are 
related to their a — > limit as follows: 

^ - ( 1 - y (^7m^m - (1 - ) ^ (24) 

where rn^ is the bare quark mass, parameter of the Wilson propagator, and is the 
lattice covariant derivative and z arbitrary parameter, corresponding to the possibility of 
using the equations of motion [^. We choose z = 

However, the systematic deviations of the order 0{g'^"'a) are generated by the action 
Sp- With the hope to diminish them we would like to eliminate tadpole diagrams within 

18 



the mean-field approximation [[l^]. The tadpole diagrams make average links differ from 
one: < f/^ >^ 1. The suggestion of fl^ is to absorb these factors by the nonperturbative 
renormalization of the parameters of the lattice action. Identifying < > = < Un >^^'^ 
in the mean-field approximation, we devide every hnk in (|2|) by < Ua >^'^- This 
implies that we use a new k = k\y / < Un >^^^ in Sp and nontrivial value of the coefficient 
c against clover term (|2^) : c = < Un >~^^^. 



2.3 Blocked lattice 

We have obtained a set of blocked SU (3) configurations from the QCD-TARO collabo- 
ration [jl5|. Starting from a 32^ x 64 lattice at /3 = 6, corresponding to ~ 2GeV 
1%, two blocking steps were performed. This brought the lattice size to 8^ x 16, with a 



rather large lattice spacing ~ 500MeV. The blocking process preserves information 
about large distances, but of course does not allow us to study short-distance effects. In 
| T8[| the clover- ||19|-and-tadpole-improved |jl6[ quark propagators were computed on these 



blocked configurations. These propagators were used in to evaluate the moments of 
the correlator of vector currents. Here we fit those moments by the moments of a single 
resonance computed in the same box in accordance with section 1.1. 

Fig. 8 shows the moments of the correlator of pseudoscalar currents. Three sets 
of data correspond to three different values of n = {0.1111,0.1031,0.0983}. Statistical 
errors are shown on the graph. The upper curve, corresponding to relatively light quarks 
K = 0.1111, {rripion ~ SOOMel^) comes from 28 configurations. The other two sets of 
data correspond to heavy quarks, where statistical errors become small for already ~ 10 
configurations. The solid lines connect the ratios of the moments of a single resonance, 
computed on the same lattice size 8'^ x 16. The input mass of that resonance m^es is a 
fit parameter, which we find to be: rrires = {1-08, 1.95, 2.36}. In this approach, the fitted 
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resonance mass automatically corrects for the systematic deviation (^Ol) . One can see that 
the solid lines fit the Monte-Carlo data perfectly at all points despite the nontrivial shape 
of the curves. (The first two ratios ri^2 are not to be considered because they are subject 
to subtractions as explained in the previous section.) This implies that the Monte-Carlo 
data on the correlator is saturated entirely by a single resonance, which is of no surprise: 
the blocked lattice sees large distances only. 

2.4 Fine lattice 

In order to be sensitive to the high-energy part of spectrum, which controls the behavior 
of the correlator at short distances, we have generated SU (3) Monte-Carlo configurations 
on a small 8^ * 16-lattice at /3 = 6 (a^^ ^ 2GeV). Quark propagators were computed for 
the clover-and-tadpole-improved action. 

Fig. 9 shows Monte-Carlo data from 20 configurations for the moments of the pseu- 
doscalar (squares), vector (octagons) and scalar (diamonds) correlators for k = 0.1000. 
Statistical errors (shown in Fig. 9, 10, 11) are rather small because the quarks are heavy. 
3000 sweeps were used for thermolization and 1000 sweeps separated configurations for 
decorrelation. No correlations from one configuration to the next were noticed. 

We fit the higher moments n > 8 by the moments of a single resonance, computed in 
an 8^ X 16 box (solid lines). Fig.9 shows three sets of double solid lines, corresponding 
to the three channels: pseudoscalar, vector and scalar. The lines are doubled to show 
the impact of statistical errors, estimated by jack-knife, on the determination of the 
resonance mass. The masses of the resonances in units of the lattice spacing are shown 
in the tables. The spectrum of the charmonium ground states thus obtained is shown in 
Fig. 10. It is gratifying to reproduce the charmonium spectrum with such a reasonable 
accuracy on a small lattice 8'^* 16. Current simulations of heavy systems are performed 
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on significantly larger lattices. The reason for this success is that the charm-anticharm 
bound states are small objects: the size of the J/tp-mesow is ~ 600Me\^. Also the 
splitting between S-wave states and P-wave states has a perturbative nature within the 
nonrelativistic approximation, and is accounted for by the leading loop of free quarks. 

The new feature that one observes in Figs. 9, 11 is that the lower moments deviate 
from the single-resonance curve. This was not seen at all on the blocked lattice. We 
attribute this fact to the short- distance contributions. 

We have modelled the hadronic spectral density as in with the phenomenologi- 
cal parameters of the charmonium spectral density ( |TTD and computed the moments of 
the corresponding correlator on an 8^ x 16 lattice in accordance with (|l^). The phe- 
nomenologically obtained lattice moments, shown in Fig. 11 by the dashed line, fit the 
Monte-Carlo data perfectly. It is obviously the contribution of the high-energy part of 
the hadronic spectrum which accounts for the deviation of the Monte-Carlo data from the 
single-resonance curve at lower moments. The contribution of the continuum spectrum 
in the smaller ratios is obviously quite large, as it accounts for the strong deviation of the 
Monte-Carlo data from the single-resonance curve. In contrast to that, the continuum 
contribution is quite small for n > 6 as compared to the resonance one. The helps to 
determine the resonance mass with better accuracy. 

Fig. 12 compares the two ways of extracting the resonance mass from Monte-Carlo 
data. The upper curve shows fit to the moment ratios with the single-resonance ap- 
poximation to the correlator. The curve flattens at higher moments corresponding to 
large distances, where the single resonance dominates. The conventional fit of the zero- 
momentum component of the correlator to cosh[{t — L/2)m^es\ is shown on Fig. 13. We do 
not see a wide horizontal plateau at large n on Fig. 12 (or at large t on Fig. 13) because 
the volume of the box 8^ * 16 is small and the contribution of the high-energy part of the 



21 



spectrum is noticeable. The bottom curve on Fig. 12 takes into account the contribution 
of the high-energy part of the spectrum. We use ansatz in the dispersion relation ([T3| ) 
with fixed residue of the resonance and the ratio of the effective continuum threshold to 
the resonance mass. Those two parameters are taken to have phenomenological values as 
in ([TT|). The resonance mass then remains a fit parameter. The procedure is applied to the 
pseudoscalar and scalar channels with similar values of the phenomenological dimension- 
less spectrum parameters / and So/m^g^. Our typical statistical error in the vector and 
pseudoscalar resonance masses is ~ ±1%, while in the scalar channel it is ~ ±5%. The 
values of the resonance mass obtained this way are shown in the tables. The continuum 
contribution helps to stabilize the resonance mass at smaller n. On bigger lattices one 
can probe longer distances, where the continuum contribution is very small, and obtain 
stable curves for the resonance mass in a wide range of n without using the continuum 
contribution. 

The model of the spectrum at high energies in the form of a smooth perturbative 
continuum is, in principle, better than the two-resonances fit because this model respects 
the global duality, which requires that the spectral density does not decrease at high en- 
ergies. In the world with no light quarks (corresponding to the quenched simulations) the 
smooth continuum models an infinite (in the limit of infinite cutoff) series of resonances, 
while in the world with light quarks the effective continuum threshold is close to the 
D-D-meson pair production in the charmonium vector channel. Since the model of the 
smooth perturbative continuum is somewhat rough in any case, one should not consider 
bad quantities, which are very sensitive to the precise form of the continuum spectrum. 
Good quantities are those related to intergrals over the spectral density of the form of 
eqn.(ll). 

The value k = 0.1060 reproduces m,,^ = 2.979GeV, for the lattice spacing a'^ = 
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l.QOGeV^, and it reproduces mj/^ = S.OOTGeV for a ^ = 1.92GeV. These values of 



the lattice spacing are consistent with the string-tension measurements at /5 = 6 [p!7[| . 
Neighboring values of k are also possible, if one changes the lattice spacing accordingly; 
to leading order the changes in k and in a compensate each other, and our physical results 
on the charmonium spectrum and the charmed quark mass are not affected. The value 
of the mass splitting within the S'-wave that we obtain with the tadpole improved clover 
action is mj/^ — m^^ = 87 ± 9 MeV. 

Fig. 11 is analogous to Fig. 1. It demonstrates the important fact that the high- 
energy part of the spectrum and the low-energy modes can be seen on the same lattice 
of limited size. In the next section we study the lower moments. We fit them with the 
corresponding expressions of perturbative QCD. 



3 The renormalized heavy-quark mass 
3.1 Off-shell mass on the lattice 

While the large- distance behavior (n ^ 1) of the correlator (0) is dominated by a single 
resonance, the short- distance behavior (n ~ 1) is determined by one-loop of heavy quarks. 
Since the Monte-Carlo simulation is done in the lattice theory with noticeable finite- 
volume and cutoff effects, we have to fit Monte-Carlo data with the loop-expansion of 
perturbative QCD evaluated in the lattice theory as well. In other words, we have to 
know the systematic deviations in the coefficients a„(a, L), 6„(a,L), c„(a,L) of eqn.(|TUp. 
As the two-loop cts-correction is small in r2,3,4 we first concentrate on one loop of free 
quarks and clarify its systematic deviations completely. 
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Consider the propagator of free Wilson fermions: 



S{q) 



m" 



Ylr {I- cos{q^)) + ml 



(25) 



At the tree level the fermion mass is related to the hopping parameter k in a simple 
way: 



am. 



1 

2k 



- 4r 



(26) 



The parameter r is in the interval < r < 1. We use the notation to emphasize that 
this mass is a parameter of the lattice theory. 

Evaluating numerically the moments of the correlator (^ in the free-fermions ap- 



proximation (with Wilson propagators (|25|) for r = 1), we obtain the ratios r„ as shown 
on Fig. 14 for periodic boundary conditions. Again, as in the case of the resonance, 
one observes a strong volume dependence: higher moments, corresponding to long dis- 
tances, feel the boundary of the box. As the volume increases S'* 32^ the function 
Tn approaches the form of the continuum theory shown in Fig. 1 by diamonds. In the 
continuum theory the moment ratios of vector, pseudoscalar and scalar correlators in the 



one-loop approximation are known to be |T0 



n 



^psc 



m? + ?>n/2 Ami 
n-l 1 



n + 1/2 Ami 
n-l 1 



ra + 3/2 Ami 



with the universal property: 



Am'^ 



(27) 
(28) 
(29) 

(30) 



where mc is the tree-level quark mass in the continuum theory. On Fig. 15 we compare 
the numerically computed moments of the pseudoscalar correlator for free Wilson quarks 



on the 32*^ lattice for k = 0.1, 



1) (solid line) with the corresponding values of the 
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continuum theory (^) (diamonds) for rhf. = 1. The discrepancy between the sohd hne 
and the diamonds is too large to remain unexplained. The asymptotics of the solid line 
at large n, determined by the large-distance {q 0) behavior of the correlator, deviates 
from (^0|). To understand this we explore the fermion propagator (^) at small momenta. 
In the limit q ^ one finds the following behavior of the denominator in (^51): 



1 2 



J2r (1 -cos(g^)) + m"^ 



q' {1 + r-m^) + {m^Y (31) 



w 

C 5 



One can see here the lattice renormalization of the fermion wave-function ^/l + ram] 
reported previously |jT6|, as well as the following tree- level lattice renormalization of the 



fermion mass: 

1 1 



■n^oo • - X2 At w\2 + ^ tt m"^ ) (32) 

Eqn. (|3^ ) is 0{a) systematic deviation. The fermion mass of the continuum theory is 
related to k in a way more complicated than (p6D even at the tree level. The ratios r„ of 
the continuum theory corresponding to mdm^ = 1) are shown in Fig. 15 by the dashed 
line. Deviation of the solid line from the dashed line at n ~ 10 is a finite- size effect, which 
would disappear on lattices larger than 32^. The equality of the asymptotics (n — > oo) 
of the solid line and the dashed line is expected by construction. The remarkable overlap 
between the solid and dashed lines over a wide range of n is less trivial. It means that 
the correlator (^ composed of free Wilson fermions behaves as in the continuum theory 
over a wide range of distances, provided the relation (^) between the masses rhc and 
is used. Fig. 16 illustrates how well this works for different k. 

The tree-level relation (|3^ ) is exact to all orders in a". This is because the mass rhc 
is the off-shell mass, defined at = in contrast to the pole mass. As promised in the 
introduction, we never use the pole mass. 

Figs. 14 and 16 show the moment ratios of the subtracted correlator as well. The 
corresponding lines terminate at n = 2. Again, just like in the case of a single resonance 
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(Figs. 2,3,4) the second ratio r2 is affected by the subtraction. The effect is seen to be 
smaller at smaller values of the fermion mass (higher n). The third ratio r-^ is affected 
very little, while the moment ratios for n > 3 are insensitive to subtractions. 



3.2 Fit of the lower moments 

As explained in the section 1, the correlator of heavy-quark currents is special in the way 
that the lower moments, which originate from the short- distance part of that correlator, 
are determined by perturbative QCD. The lowest moments must be saturated by one-loop 
of free quarks. Therefore we fit those moments, obtained by Monte-Carlo on the lattice, 
with one loop of free Wilson quarks. 

In this section we neglect 0{as) terms and obtain an approximate value of the renor- 
malized heavy-quark mass rhc from the ratio r^. 



« \^ (33) 



We use the coefficient of the continuum theory because the coefficients a„ of the con- 
tinuum theory (pTf) - (El) fit the lattice correlator, composed of Wilson propagators, so 



well (see Fig. 15, 16), provided the mass relation (|3^) is used. 

Fig. 18 shows the fit of Monte-Carlo data at n = 3 by the moments of the pseudoscalar 
correlator, saturated by one loop of effectively free Wilson quarks. Those quarks are found 
to have k = 0.0994, which for free Wilson quarks implies, am^ = 1.02 because of eq. (^6]) . 
From the previous section we know that for such a large mass it is important to correct 



for the 0{a) systematic deviation ([321) . One obtains amc{am'^ = 1.02) = 0.72, which for 
= 1.9 GeV gives: 

rhc = 1-37 ± 0.01 GeV (34) 
where the error is associated with the contribution of the high-energy part of the spectrum. 
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Other channels, such as vector or scalar, are just as good as the pseudoscalar one 
for the determination of the renormalized quark mass. The ratio r4 would be also just 
as good as to determine rhc, since the relation between and r4 is calculable within 
perturbative QCD. The difference is O ( a s) -effect. We consider it in section (3.4). 

Note, that the relation (|32|) is derived in the infinite volume. Since we study rather 
small lattice 8^ * 16 we would like to be more accurate with the finite- volume effects. The 
moments of the continuum theory ( p8|) agree with the lattice moments (see Fig. 15, 16), 
computed in a rather large box 32^. In the mean time one can notice in Fig. 14 that the 
moments of the small lattices like 8'^ * 16 deviate from the asymptotic L —>■ oo behavior 
even at small n = {1, 2, 3}. In the next section we take into account this deviation. 

3.3 The lattice dispersion relation for the fermionic continuum 

We have seen on Fig. 15, 16 that the moment ratios r„ of the lattice correlator in the 
free- Wilson-quarks approximation follow the corresponding ratios of the continuum theory 
with the correct value of the quark mass. The moments ( pT]) - ( pOf ) of the continuum theory 
can be calculated with the help of the spectral represenation (|^) with the spectral density 
corresponding to one-loop of free quarks. For the three channels under discussion one can 
use: 

p-^(.) = ^ {s + 2{m,r) ^1 - 4(m,)V. (35) 
//-(.) = A , ^1 _ 4m,/. (36) 
p-{s) = {s- 4(mj2) - 4(m,)V. (37) 

We would like to see how well the lattice variant of the dispersion relation ([131) works 
for the continuum spectrum of free fermions. We can use the spectral densities (^) - 
(p^) in the eqn. (|13[) and compare the result with the correlator composed of free Wilson 
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propagators. We plot the moment ratios obtained in those two ways on Fig. 17 in the 
case of small lattice 8^ * 16. We fit the moments of the lattice correlator of free Wilson 
fermions with corresponding moments of the correlator computed via the dispersion rela- 



tion (13) with the spectral density (p6|). This is a one-parameter fit. That parameter is 
the quark mass of the continuum theory m^. We find the agreement illustrated by Fig. 
17 remarkable. 

Thus, the use of the lattice dispersion relation (0) for the continuum spectrum of 
quark- antiquark pair helps to fit the moments on small lattices. It appeares as an efficient 
instrument to take into account-finite size effects. On bigger lattices one can simply use 
the formula fB2l). 



3.4 Fixing the subtraction scheme 

In this section we estimate the value m^^'^(mc) of the renormalized charmed-quark mass 
in the minimal subtraction scheme normalized at the Euclidean point fi = rrtc. In order 
to fix the subtraction scheme one should be sensitive to the a^-correction in the ratios 
r„. The study of the systematic deviations of the coefficient 6„ is beyond the framework 
of the present paper. However, since the term 6„ as is small in the lower ratios ^"2,3,4 we 
will take its values from the continuum theory. This is our approximation. To be more 
precise, we write the lower ratios as 

Or,. . 



Ami 

i = 1 + ashn/an (38) 

We take the value of the strong coupling constant asirric) ~ 0.3 |. Based on the two- loop 
calculations of the continuum theory |]10[ we obtain for the vector, pseudoscalar and scalar 



^To fix precisely the normalization point of as in every moment one needs to know the three-loop 
correction in the ratios (pO|), which is not known at present 
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channels: 

er = 1-02 , a^'^ = 0.95 , er" = 0.97 

Cf" = 0.96, Cf^ = 0.92, ^f« = 0.91 (39) 

The a^-term is seen to introduce a few-percent correction to the ratios r2,3,4. Therefore, 
the presently unknown cutoff effects in the coefficient 6„ have a small effect on the value 
of the quark mass fhc- 

We fit the lower moment ratios rs^4 of our Monte-Carlo data using the lattice spectral 
representation, as explained in the previous section, to relate more accurately the quark 
mass of the continuum theory fhc to the corresponding Wilson propagator in the finite- 
volume lattice theory. 

The results for the renormalized charmed-quark mass are shown in the tables. The 
tables demonstrate that the dimensionless ratio of the renormalized quark mass to the 
resonance mass fhc/mres is very stable with respect to significant variation of k. This 
implies that the renormalized heavy-quark mass is less sensitive to the uncertainty in the 
value of the lattice spacing than one might expect. 

The relation rric between the renormalized heavy-quark mass and bare param- 
eters of the lattice theory is, in principle, rather non-trivial. On Fig. 18 we plot the 
moments of one loop of the free Wilson quarks with the input k — 0.1060 (dashed fine) 
to show how different they are from the moments of one loop of free quarks with k corre- 
sponding to the renormalized mass (bottom line) . At large values of (5 the relation between 
K, and h is perturbative. We have generated 20 gauge configurations for /3 = 20 with the 
same Monte-Carlo parameters as for /3 = 6. The corresponding moments of the correla- 
tor with the clover-and-tadpole-improved propagators for the same input k = 0.1060 are 
shown on Fig. 18 by squares. The Monte-Carlo data for = 20 is fitted perfectly at all 
n by one loop of free-quarks (upper solid line) with k' — 0.1034. We find numerically the 
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relation : 

k' = K < Un{f3 = 20) >^/^ (40) 

Obviously the data at /3 = 20 shows the correlator at distances so small that the effects 
due to interaction are not seen. The contribution of the clover term is strongly suppressed. 
The relation (^) is due to the fact that we used the mean-field-tadpole-improved action 



T6[|. Generating the correlator by Monte-Carlo at various high values of jS one can study 



the correlator in the weak-coupling expansion ||I4 . 



4 Conclusions 

We have formulated in details the way to calculate the renormalized quark mass us- 
ing Monte-Carlo data on the correlator of heavy-quark currents. We separate the long- 
distance and the short- distance parts of the correlator by studying the higher and the 
lower moments of the correlator respectively. The long-distance part is used to fix the 
scale of the lattice theory, as it is saturated by the lowlying resonance. The short- distance 
part is fitted by the loop expansion of perturbative QCD. The coefficients of that expan- 
sion are contaminated by the cutoff effects and the finite-volume effects. Those system- 
atic deviations can be studied order by order in the loop-expansion. We have worked 
out systematic deviations of the one loop of Wilson quarks, which saturates the very 
short-distance part of the correlator (the lowest moments). Then fit to Monte-Carlo data 
produces the value of the renormalized heavy-quark mass. We find for the charmed quark 
mass: m^^^{mc) = 1.22{2)GeV - from the vector channel and m^^^{mc) = 1.18(2)GeV^- 
from the pseudoscalar channel. There are two sources of uncertainties in these values. 
The first is the ansence of scaling on the present lattice. The data for higher /3 > 6 as well 
as improved actions should helpful in this respect. The second uncertainty is in fixing the 
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scale. One needs bigger volume to determine the resonance mass with better accuracy. 
The dimensionless ratio of the renormalized heavy-quark mass to the resonance mass is 
stable with respect to variation of n provided the value of n corresponds to heavy quarks. 

Our central value for the charmed quark mass in the vector channel can be compared 
to the corresponding estimate of the continuum theory : m^^^{mc) = 1.236*6^0. There 
are no estimates of the charmed quark mass from the pseudoscalar or scalar correlators 
in the continuum theory, since the corresponding spectral densities are not observable 
experimentally. 

A comprehensive study of the correlator of heavy-quark currents at both large and 
small distances is a challenging problem. One principal question we tried to clarify here is 
whether one can observe simultaneously, on the lattice sizes presently available to Monte- 
Carlo simulations, both short-distance phenomena, described by asymptotic freedom, and 
long-distance tails of hadronic correlators, saturated by low-lying resonances. In order to 
see short-distance phenomena one has to keep the lattice spacing small. Then one needs 
many lattice sites in order to reach large physical scales, where the low-energy part of 
the hadronic spectrum dominates. To solve this problem it is crucial to combine the 
extension of dispersion relations to the lattice with the improved actions. The former 
helps to fit correlators in smaller volumes, while the latter helps to see QCD perturbation 



theory on coarse lattices. Further study of the a^-corrections [Tj| to the correlator of the 
heavy-quark currents on the lattice is necessary to specify the renormalization scheme of 
the extracted quark mass more accurately and measure the renormalized strong coupling 
constant. 
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Table 1. Charm-quark mass as a result of the fit of the ratios r^, of the subtracted 
correlator for k — 0.1000. The masses rrires and rhc are in units of the lattice spacing. The 
masses m^^[aji^] and m^'^[a^^] are in GeV . They are obtained assuming that the lattice 
spacing is fixed from the vector channel and from the pseudoscalar channel respectively. 
The indicated errors are statistical. 
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1.87(10) 
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Table 2. The same as in Table 1 for k = 0.1060. 
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Figure 1: Ratios of the neighboring moments in the continuum theory for 
the correlator of interpolating currents of the J / t/j-meson. Burst symbols 
correspond to the evaluation of the correlator via dispersion relation with the 
phenomenological hadronic spectrum (11). Squares show the contribution of 
a single J/r/j-resonance. Diamonds are moments of one loop of free quarks 
with mass nic = 1.26 GeV. 




2 4 6 8 10 12 

n 

Figure 2: Ratios of the neighboring moments of the single resonance ap- 
proximation to the correlator, on three lattices of different sizes: S"*, f G"*, 32'*. 
The resonance mass (m^es = 0.5) is given in the units of the inverse lattice 
spacing. Solid lines are to guide the eye. The lines becomes double near 
n = 2: then one line corresponds to the subtracted correlator (terminating 
at n = 2), the other - to the unsubtracted one. 
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Figure 6: Systematic deviation 0{{amresY) of the resonance mass m^es, ex- 
tracted from the ratios function of the input resonance mass m^es- 
The sohd hne is the fit: 1 + 0.04m^gg. 
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Figure 7: Systematic deviation (9((amres)^) of the resonance mass, extracted 
in two different ways: the dotted hne is a conventional fit of the zero- 
momentum correlator to cos/i[(t — iv/2)mres]; the dashed line is the resonance 
mass obtained from the ratios r„ . 
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Figure 8: Ratios r„ of the Monte-Carlo data for the pseudoscalar currents 
on the lattice 8"^ X 16, blocked from 32"^ X 64, /3 = 6. Solid lines show the 
one-parameter {nires) fit by the ratios of a single resonance. Three sets of 
data correspond to k = {0.1111,0.1031,0.0983}. 
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Figure 9: Monte-Carlo data on the 8"^ X 16, /3 = 6 - lattice for pseudoscalar 
(squares), vector (octogons) and scalar (diamonds) currents, n = 0.1000. 
Solid double lines show the ratios of a single resonance, whose mass is fitted 
to the higher ratios. Double lines refiect statistical errors. 
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Figure 10: Charmonium spectrum obtained on the 8"^ X 16-lattice at /3 = 6 
with the clover-and-tapole-improved action. The vector mass mj/^ = 
3.097GeV is used for normahzation. The wide horizontal lines are exper- 
imental points. 




Figure 11: Monte-Carlo data on the 8"^ X 16, /3 = 6 - lattice for vector 
currents, n = 0.1060. The solid line is a single resonance approximation. 
The dashed line is the correlator with the phenomenological charmonium 
spectrum (11) incorporated via the dispersion relation (13). 
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Figure 12: Resonance mass from the fitting the moment ratios of the Monte- 
Carlo data on the 8"^ X 16, /3 = 6 - lattice for vector currents, n = 0.1060. 
The upper line is a single resonance approximation. The bottom line cor- 
responds to the correlator with the spectrum of eqn. (5) incorporated via 
the dispersion relation (13). The residue of the resonance / and the ratio 
■^o/'m-les kept fixed and taken to have phenomenological values (11) while 
the resonance mass is a fit parameter. 
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Figure 13: Resonance mass from the conventional fit of the sero-momentum 
component of the correlator of vector currents (the same Monte-Carlo data 
as in Fig. 12) to cosh[[t — L/2)mres]- 
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Figure 14: Ratios of the neighboring moments of the subtracted (dashed 
hues) and unsubtracted (sohd hues) correlator of pseudoscalar currents, sat- 
urated by one loop of free Wilson fermions with n = 0.1 [am'^ = 1) and 
periodic boundary conditions, computed numerically on lattices of the indi- 
cated size. Lines connect the data to guide the eye. 




Figure 15: Ratios r„ for one loop of free Wilson fermions (pseudoscalar cur- 
rents), computed numerically on the lattice 32'* (solid line) for mjf [k, = 0.1). 
Diamonds are the ratios of the continuum theory for mjf = 1. The dashed 
line shows the ratios of the continuum theory with the corrected mass rfic- 




Figure 16: Same as in Fig. 9, for n = {0.1111, 0.1053, 0.1, 0.0952}, correspond- 
ing to = {0.5,0.75,1,1.25}. 




Figure 17: Fit of the pseudoscalar correlator of free Wilson quarks (solid line) 
with the dispersion relation (dashed line) for n = {0.1111,0.1060,0.1000}, 
corresponding to amjf = {0.5,0.72,1}. 
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Figure 18: Monte-Carlo data lor the pseudoscalar currents on the 8"^ X 16- 
lattice at /3 = 6 (diamonds), /3 = 20 (squares); n = 0.1060. The dashed 
line shows the ratios lor one loop ol Iree Wilson lermions with n = 0.1060. 
The solid lines show the ratios lor one loop ol Iree Wilson lermions: the 
bottom solid line is lor k = 0.0994, corresponding to the physical quark 
mass rhc = LSTGeV, the top one is lor one loop ol Wilson lermions with 
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